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Abstract 

We have searched for the lepton- flavor- violating decays r~ — > £~KS and £~K§Kg (£ = e or /i), 
using a data sample of 671 fb -1 collected with the Belle detector at the KEKB asymmetric-energy 
e + e~ collider. No evidence for a signal was found in any of the decay modes, and we set the following 
upper limits for the branching fractions: £>(r~ — >■ e~K§) < 2.6xl0~ 8 , B(r~ — > fi~Kg) < 2.3xl0 -8 , 
B(t- -»• e'KlKl) < 7.1 x 10" 8 and B(r~ -»• ^K^Kg) < 8.0 x 10~ 8 at the 90% confidence level. 

PACS numbers: 11.30.Fs; 13.35.Dx; 14.60.Fg 
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INTRODUCTION 



Lepton flavor violation (LFV) in charged lepton decays is forbidden in the Standard 
Model (SM) or highly suppressed if neutrino mixing is included. However, LFV appears 
in various extensions of the SM. In particular, the lepton-flavor-violating decays t~ — > 
£~K§ and r~ — > £~K^Kg (where £ = e or /i ) are enhanced in supersymmetric and many 
other models [H-H]. Some of these models predict branching fractions which, for certain 
combinations of model parameters, can be as high as 10~ 7 ; this level is already accessible in 
high-statistics 5-factory experiments. Previously, we obtained 90% confidence level (C.L.) 
upper limits for the t~ — > £~K® branching fractions (£>) using 281 fb _1 of data; the results 
were B[t~ ->■ e~K^) < 5.6 x 10" 8 and B(r~ -> fi'Kg) < 4.9 x 10" 8 Q. The BaBar 
collaboration has recently obtained 90% C.L. upper limits of B[r~ — > e~K®) < 3.3 x 10~ 8 
and £>(r~ — > fi~K§) < 4.0 x 10~ 8 using a data sample of 469 fb _1 (8[. The most restrictive 
existing upper limits B{r~ ->■ e~K^) < 2.2 x 10~ 6 and B{r~ -> fi~K^lQ) < 3.4 x 10~ 6 
at the 90% C.L. were set by the CLEO experiment using 13.9 fb _1 of data In this 
paper, we present a search for the lepton-flavor-violating decays t~ — ^ £ ^-s an d ^~^s^-s 
{£ = e or ji) [|] using 671 fb _1 of data collected at the T(4S') resonance and 60 MeV below 



with the Belle detector at the KEKB asymmetric-energy e + e~ collider [10 . 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel- like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL), all located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect mesons and to identify muons 



(KLM). The detector is described in detail elsewhere [11 . 

Leptons are identified using likelihood ratios calculated from the responses of vari- 
ous detector subsystems. For electron identification, the likelihood ratio is defined as 
V(e) = C e /(C e + C x ), where C e and C x are the likelihoods for electron and non-electron 
hypotheses, respectively, determined using the ratio of the energy deposit in the ECL to the 
momentum measured in the SVD and CDC, the shower shape in the ECL, the matching 
between the position of the charged track trajectory and the cluster position in the ECL, 
the hit information from the ACC, and the dE/dx information in the CDC [l2|. For muon 
identification, the likelihood ratio is defined as V{\i) = C h J{C^ + C n + Ck), where £ At , C w 
and Ck are the likelihoods for the muon, pion and kaon hypotheses, respectively, based on 
the matching quality and penetration depth of associated hits in the KLM [l3|. For this 
measurement, we use hadron identification likelihood variables based on the hit information 
from the ACC, the dE/dx information in the CDC, and the particle time-of-flight from the 
TOF. To distinguish hadron species, we use likelihood ratios, V(i/j) = Ci/(Ci + where 
Li (£j) is the likelihood for the detector response to a track with flavor hypothesis i (j). 

In order to optimize the event selection and estimate the signal efficiency, we use Monte 
Carlo (MC) samples. The signal and background events from generic t + t~ decays are gener- 
ated by KKMC/TAUOLA [14| . The signal MC samples are generated by KKMC assuming 
a phase space model for the r decay. Other backgrounds, including BB and continuum 
e + e~ — > qq (q = u, d, s, c) events, Bhabha events, and two-photon processes, are generated 



by EvtGen [151 ]. BHLUMI |16|, and AAFH [17J, respectively. The Belle detector response 



[|] Unless otherwise stated, charge-conjugate decays are included throughout this paper. 
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is simulated by a GEANT 3 18[ based program. The event selection is optimized mode- 
by-mode since the backgrounds are mode dependent. All kinematic variables are calculated 
in the laboratory frame unless otherwise specified. In particular, variables calculated in the 
e + e~ center-of-mass (CM) system are indicated by the superscript "CM" . 



DATA ANALYSIS 

We search for t + t~ events, in which one r (signal side) decays into £K$ or £K$Kg, 
while the other r (tag side) decays into a final state with one charged track, any number of 
additional photons and neutrinos. We reconstruct each K$ meson candidate from a 7t + tt~ 
pair. By selecting decays into one charged track on the tag side, we reduce background 
from BB and qq events. All charged tracks and photons are required to be reconstructed 
within a fiducial volume, defined by —0.866 < cos 9 < 0.956, where 9 is the polar angle with 
respect to the direction opposite to the e + beam. We select charged tracks with momenta 
transverse to the e + beam p t > 0.1 GeV/c and photons with energies > 0.1 GeV. 

Candidate r-pair events are required to have four or six charged tracks with zero net 
charge for the £K§ and £K^K^ modes, respectively. Events are separated into two hemi- 
spheres corresponding to the signal (three-prong and five-prong for the £K§ and £K^K^ 
modes, respectively) and tag (one-prong) sides by the plane perpendicular to the thrust 
axis [l9j. 

We require one or two Kg candidates for the £K§ and £K§K§ modes, respectively. The Kg 
is reconstructed from two oppositely charged tracks on the signal side that have an invariant 
mass 0.482 GeV/c 2 < M^+ n - < 0.514 GeV/c 2 , assuming the pion mass for both tracks. The 
7r + 7r~ vertex is required to be displaced from the interaction point (IP) in the direction of the 
pion pair momentum [2(j. In order to avoid fake K$ candidates from photon conversions 
(i.e., 7 — > e + e~), the invariant mass reconstructed by assigning the electron mass to the 
tracks, is required to be greater than 0.2 GeV/c 2 . The electron and muon identification 
criteria are V(e) > 0.9 with momentum p > 0.3 GeV/c and V(fi) > 0.9 with p > 0.6 GeV/c, 
respectively. In order to take into account the emission of bremsstrahlung photons from the 
electron, the momentum of each electron track is reconstructed by adding the momentum 
of every photon within 0.05 rad of the track direction. The electron (muon) identification 
efficiency for the £K% modes is 92% (87%) and that for the £KgK% modes is 79% (81%). The 
difference of efficiencies between £K$ and £KgK% is due to the different signal momentum 
distributions. The probability to misidentify a pion as an electron and a muon is below 0.5% 
and 3%, respectively. 

In order to suppress background from qq events, the following requirements on the number 
of the photon candidates on the signal and tag side (n^ IG and n^ AG ) are imposed: n^ 1G < 1 
and n^ AG < 3. For the £Kg modes only, we also require n^ AG < 1 if the track of the tag 
side is a lepton to reduce the background, in particular from D + — > i + uKg(— > 7r°7r°). 

To ensure that the missing particles are neutrinos rather than photons or charged particles 
that fall outside the detector acceptance, we impose additional requirements on the missing 
momentum vector, p m i SS , calculated by subtracting the vector sum of the momenta of all 
tracks and photons from the sum of the e + and e~ beam momenta. We require that the 
magnitude of p m i S s be greater than 0.4 GeV/c and that its direction point into the fiducial 
volume of the detector. Since neutrinos are emitted only on the tag side, the direction of 
p m i SS should lie within the tag side of the event. The cosine of the opening angle between 
p m i SS and the tag-side track in the CM system, cos^^L^g, should be 0.0 < cos0f a g_ miss 
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FIG. 1: Kinematic distributions used in the event selection of the r~ — > fi~Kg mode: (a) the 
cosine of the opening angle between a charged track on the tag side and the missing momentum 
in the CM system (cos ^tagf-niiss) > an< ^ 3 ) the magnitude of the thrust. The signal MC (r~ — > 
fi~K§) distributions with arbitrary normalization are shown for comparison; the background MC 
distributions are normalized to the data luminosity. Selected regions are indicated by the arrows 
from the marked cut boundaries. 
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FIG. 2: Kinematic distributions used in the event selection of the t~ — > ^K^K^ mode: (a) the 
cosine of the opening angle between a charged track on the tag side and the missing momentum 
in the CM system (cos ^gf-niiss) i anc ^ 3 ) the magnitude of the thrust. The signal MC (r~ — > 
/i~KgKg) distributions with arbitrary normalization are shown for comparison; the background 
MC distributions are normalized to the data luminosity. Selected regions are indicated by the 
arrows from the marked cut boundaries. 

for both modes (see Figs. [JJ (a) and [2] (a)). For the £K$ modes, we also require that 
cos #° a g_ miss < 0.99 to reduce background from Bhabha, and two-photon events, as 

radiated photons from the tag-side track result in missing momentum if they overlap with 
the ECL clusters associated with the tag-side track. 

To reject the qq background, the magnitude of the thrust is required to be larger than 
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TABLE I: The selection criteria for the missing momentum (p m i ss ) and missing mass squared 
( m miss) correlations, where p m i ss is in GeV/c and m ^ i sR is in (GeV/c 2 ) 2 . 



Modes 


Hadronic tag 


Leptonic tag 


/Kg 


Pmiss > -3.0 x - 0.9 
p miss > 3.5 x m miss - 1.1 


Pmiss > -4 x m miss - 1.0 
Pmiss > 1-8 x m miss -0.8 


£K°K° 


Pmiss > -2 x m miss - 1.0 
p miss > 2 x m miss -1.0 


Pmiss > -2 x m miss - 1.0 
Pmiss > 1-3 x m miss - 0.8 



0.9 (see Fig. [T](b) and[2] (b)). The invariant mass reconstructed from the charged track and 
any photon on the tag side is required to be less than 1.0 and 1.777 GeV/c 2 for the £K§ and 
£K®K® modes, respectively. For the £K§ modes, we impose a kaon veto £(K/tt) < 0.6 if the 
track on the tag side is a hadron, to suppress e + e~ — > qq background; due to the conservation 
of strangeness by the strong interaction, the Kg in such events is often accompanied by 
another kaon. 

All kinematic distributions for the £Kg modes shown in Fig. [1] are in reasonable agreement 
between data and background MC while those for the iKgK^ modes shown in Fig. [2] clearly 
differ. This difference between the data and background MC in Fig. [2] originates from our 
poor knowledge of the branching fractions B{t~ — > 7r~KgKgZ/ T ) = (2.4 ± 0.5) x 10~ 4 and 
B(t~ -> n-K°K°ii v T ) x (B(K° -»■ Kg)) 2 = (3.1 ± 2.3) x 10~ 4 x 1/4 and the dynamics 
of these decays |2l| . Since the final estimate of the background uses information from the 
data, this discrepancy does not directly affect our results. 

Finally, to suppress backgrounds from generic t + t~ and qq events, we apply a selection 
based on the magnitude of the missing momentum p m iss and the missing mass squared m miss . 
The latter is defined as K miss — p miss , where K m iss = K tot ai — K vis , K tota i is the sum of the beam 
energies and K vis is the total visible energy. We apply different selection criteria depending 
on the type of one-prong tag: the number of emitted neutrinos is two if the tagging track is 
an electron or muon (leptonic tag) while it is one if the tagging track is a hadron (hadronic 
tag). The requirements are listed in Tabled] (see also Fig. [3]). 

RESULTS 

Signal candidates are examined in two-dimensional plots of the £K§ and iK^Kg invariant 
mass, M S ig (= M eK o, M iK o K o), and the difference of their energy from the beam energy in 
the CM system, AE. A signal event should have M s ; g close to the r-lepton mass and AE 
close to zero. For both modes, the M S i g and AE resolutions are parameterized from the 
MC distributions around the peak region using asymmetric Gaussian shapes to take into 
account initial state radiation. The widths of these Gaussians are shown in Table [Til 

To evaluate the branching fractions, we use an elliptical signal region that contains 90% of 
the signal MC events satisfying all selection criteria. The shape of the signal region is chosen 
to minimize its area and therefore obtain the highest sensitivity. We blind the data in the 
signal region until all selection criteria are finalized so as not to bias our choice of selection 
criteria. Figure H] shows scatter-plots for data events and signal MC samples distributed 
over ±20cr in the M sig — AE plane. As MC simulation shows, the dominant background in 
the signal region comes from events with a fake lepton from a pion. Therefore, we estimate 
the number of expected background by multiplying the number of data events in the signal 
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FIG. 3: Scatter-plots of p m i ss vs. m^ isfi : (a), (c) and (e) show the signal MC (r~ — > fi~K^), the 
generic r + r~ MC and distributions, respectively, for the hadronic tags while (b), (d) and (f) 
show the same distributions for the leptonic tags. Selected regions are indicated by lines. 



region with selected hadrons (V(£) < 0.9) by the fake lepton ratio. The latter is calculated 
as the number of events in the data with P{£) > 0.9 divided by the number of events in the 
data with P(£) < 0.9 in the sideband region. For the £K§ modes we define the sideband 
region as the box inside the two horizontal lines (see Fig. [4] (a) and (b)) with the signal region 
excluded since real leptons from D + — > £ + i>K^ populate the region below the AE signal one. 
For the £K^K^ modes, events that lie within a ±20a region but outside the signal region 
are treated as sideband events (see Fig. [4] (c) and (d)). The final signal efficiency and the 
number of expected background events in the signal region for each mode are summarized 
in Table ED 

The dominant systematic uncertainties on the detection sensitivity come from K$ recon- 
struction and tracking efficiencies. These are 4.5% per Kg candidate and 1.0% per track. 
Other sources of systematic uncertainties are: lepton identification (2.2—2.7)%, MC statis- 
tics (0.8 — 1.0)%, trigger efficiency (0.01—0.4)%, and integrated luminosity (1.4%). The un- 
certainty from B(Kg — > 7T + 7r _ ) is negligible. All these uncertainties are added in quadrature 
to provide total systematic uncertainties that range from 6.6% to 11.3%. 

Finally, we examine the blinded region and find no data events in the signal region for 
any of the decay modes (see Fig. HJ). Therefore, we set the following upper limits on the 
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TABLE II: Summary of M sig and AE resolutions (MeV/c 2 ) and a^f (MeV)). Here 

^high ^low) means 

the standard deviation on the higher (lower) side of the peak. 



Mode 


high 


low 


high 


low 
a AE 




7.3 


7.5 


19.4 


30.0 




6.2 


6.8 


19.1 


26.4 




5.6 


6.4 


12.6 


21.9 




5.2 


6.0 


11.2 


17.2 



> 

o 
C3 



. 0-2 



-0.4 



(a) t~ 



■ e - K°. 




1.65 1.7 1.75 1.8 1.85 J A, 

M eKs (GeV/c?) 



o 



(b) t-^/x-K 



,0.2 



□ □ 



1.65 1.7 1.75 



3 , 1.85 . ,1 v £ 

;Ks (GeV/c 2 ) 




! J- 85 , *. 
M, K , K , (GeV/c 2 ) 



G 
O 



(d) t' -^/u,'K°,K° 



,0.1 



-0.2 




1.7 1.75 1.8 



(Ge'v7c 2 ) 



FIG. 4: Scatter-plots of data in the M s i g - AE plane: (a), (b), (c) and (d) correspond to the ±20cr 
area for the r~ — > e _ .Kg, r~ — > yTK§, t~ — > e~K^Kg and r~ — > fi~KgK^ modes, respectively. 
Data is indicated by the solid circles. The filled boxes show the MC signal distribution with arbi- 
trary normalization. The elliptical signal regions shown by the solid curves are used for evaluating 
the signal yield. In (a) and (b), the region between the horizontal solid lines excluding the signal 
region is used as a sideband. 
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TABLE III: The signal efficiency (e), the number of the expected background events (-/Vbg) es_ 
timated from the sideband data, the total systematic uncertainty (cr sys t), the number of observed 
events in the signal region (-/V b s ), 90% C.L. upper limit on the number of signal events includ- 
ing systematic uncertainties (sgo) and 90% C.L. upper limit on the branching fraction for each 
individual mode. 



Mode 


e(%) 


A^BG 


0" S yst (%) 


N ohs 




B (xHT 8 ) 


T~ -»• e~K% 


10.2 


0.18±0.18 


6.6 





2.25 


<2.6 


r- fi-Kl 


10.7 


0.35±0.21 


6.8 





2.10 


<2.3 


T~ -> e-K%K% 


5.82 


0.07±0.07 


11.2 





2.44 


<7.1 


r- ^K$Kl 


5.08 


0.12±0.08 


11.3 





2.40 


<8.0 



branching fractions based on the Feldman-Cousins method 22J. The 90% C.L. upper limit 
on the number of signal events (sgo) is obtained using the POLE program 23], based on 
the number of expected background events, observed data and the systematic uncertainty. 
The upper limit on the branching fraction is then given by 

Sir- -> i-KKKD) < 2EB{iq ^ rYNt , (1) 



where e is the signal efficiency, B(K§ — > n + ix~) = (69.20 ± 0.5)% [2l|, and n is 1 and 2 for 
the £K® and £K^K^ modes, respectively. The value N TT = 6.17 x 10 8 is obtained from the 
product of the integrated luminosity and the cross section of r-pair production 0.919 ±0.003 
nb {24J. The resulting upper limits on the branching fractions at the 90% C.L. are 

B(r~ -> e'Kl) < 2.6 x 10~ 8 , 
B(t~ -> iTKjl) < 2.3 x 10" 8 , 
B(t~ e'K^Ki) < 7.1 x 10" 8 , 
B{j- -> [TKlKl) < 8.0 x 10~ 8 . 

For the £K§ modes, these results improve the existing upper limits by about a factor of 
2, compared to our previously published limits 0. For the iK^R® modes, these results 
improve the upper limits by factors of 31 and 43 for the eK^R® and fiK^K^, respectively, 
compared to the previously published limits obtained by the CLEO experiment [9J. 



SUMMARY 

We have searched for the lepton-fiavor-violating decays r~ — > £~K<$ and £~K^Kg {£ = 
e or fi) using data collected with the Belle detector at the KEKB e + e~ asymmetric-energy 
collider. We find no signal for any decay modes. The following upper limits on branching 
fractions at the 90% confidence level are obtained: B(t~ — > e~K%) < 2.6 x 10~ 8 , B(t~ — > 
fi-K%) < 2.3 x 10~ 8 , B(t~ -)■ e~K%Kl) < 7.1 x 10~ 8 and B(t~ -> yrK%K%) < 8.0 x 10~ 8 . 
These results are currently the most stringent upper limits for the £K§ and the £KgK$ 
modes. These limits can be used to constrain new physics scenarios beyond the Standard 
Model. 
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